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Abstract 

A commonplace view of pressure-driven turbulence in pipes and channels is as "cas- 
cades" of streamwise momentum toward the viscous layer at the wall. We present in this 
paper an alternative picture of these flows as "inverse cascades" of spanwisc vorticity, in 
the cross-stream direction but away from the viscous sublayer. We show that there is a 
constant spatial flux of spanwise vorticity, due to vorticity conservation, and that this flux 
is necessary to produce pressure-drop and energy dissipation. The vorticity transport is 
shown to be dominated by viscous diffusion at distances closer to the wall than the peak 
Reynolds stress, well into the classical log-layer. The Perry-Chong model based on "rep- 
resentative" hairpin/horsehoe vortices predicts a single sign of the turbulent vorticity flux 
over the whole log-layer, whereas the actual flux must change sign at the location of the 
Reynolds-stress maximum. Sign-reversal may be achieved by assuming a slow power-law 
decay of the Townsend "eddy intensity function" for wall-normal distances greater than 
the hairpin length-scale. The vortex-cascade picture presented here has a close analogue in 
the theory of quantum superfluids and superconductors, the "phase slippage" of quantized 
vortex lines. Most of our results should therefore apply as well to superfluid turbulence in 
pipes and channels. We also discuss issues about drag-reduction from this perspective. 



1 Introduction 



It is quite standard to view turbulence in pipes, channels and ducts, as well as other wall- 
bounded turbulent boundary layers, as spatial "cascades" of momentum in the wall-normal or 
cross-stream direction. For example, we may quote from a well-known monograph pj: 

"There exists a close analogy between the spatial structure of turbulent boundary 
layers and the spectral structure of turbulence. At sufficiently large Reynolds num- 
bers, the overall dynamics of turbulent boundary layers is independent of viscosity 
just as the large-scale spectral dynamics of turbulence is. In the wall layer of a 
turbulent boundary layer, viscosity generates a 'sink' for momentum, much like the 
dissipative sink for kinetic energy at the small-scale end of the turbulence spec- 
trum. In particular, the asymptotic rules governing the link between the large-scale 
description and the small-scale description lead to the closely related concepts of 
an inertial subrange in the turbulence energy spectrum. ..and an inertial sublayer in 
wall-bounded shear flows." — Tennekes & Lumley (1972) 

In this point of view, the logarithmic mean velocity profile is analogous to the Kolmogorov —5/3 
power-law energy spectrum, and the constant Reynolds stress through the inertial sublayer is 
analogous to the constant energy flux through the inertial subrange. As expressed in the quote 
above, high momentum from the outer region of the boundary layer is transported inward by 
the constant Reynolds stress and ultimately transferred to the wall by viscous stress. 

There is, however, another complementary process in wall-bounded turbulence of at least 
equal importance. Vorticity generated at the wall is transported outward, diffused first by 
viscosity and subsequently advected by turbulent velocity fluctuations. Due to the fundamental 
hydrodynamic principle of vorticity conservation |2l!3], the wall-normal flux of spanwise vorticity 
is constant over the whole cross-section of the flow. The "sink" of vorticity at the outer range 
is provided by cancellation with vorticity of the reverse sign transported from the wall on the 
opposite side of the flow. There is thus a spatial "inverse cascade" of vorticity co-existing with 
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the "forward cascade" of momentum. Furthermore, it may be shown that the constant flux of 
spanwise vorticity in this "inverse cascade" is crucially related to the energy dissipation and 
wall drag in turbulent pipe and channel flows. This relation was first observed in early work of 
Taylor |i4 ; , but seems not to have been as widely appreciated as it deserves in the fluid dynamics 
community. Although there have been very many studies of the Reynolds stress — or turbulent 
momentum transport — in pipe and channel flow, there have been far fewer investigations of 
the turbulent vorticity-transport tensor in classical fluid turbulence. Notable exceptions are 
the experimental works by Klewicki and collaborators [3 El El E] and a theoretical and DNS 
study of Bernard [9]. In the meantime, Taylor's observation of the relation between cross- 
stream transport of vorticity and downstream pressure drop has been rediscovered in the field 
of quantum superfluids, where it goes by the name of the "Josephson- Anderson relation" |1U | I11]. 
The cross-stream motion or phase slippage of quantized vortex lines is widely recognized to be 
a key mechanism of energy dissipation in quantum superfluids and superconductors [121 US] • 

The purpose of this paper is to expound in detail the picture of pipe and channel flows, 
or other pressure-driven turbulent flows, as spatial "inverse cascades" of vorticity co-existing 
with the more commonly recognized "forward cascade" of momentum. Although several of the 
results presented here are to be found in the existing literature, there seems to be no systematic 
presentation of this fundamental point of view of wall-bounded turbulence. It is hoped that 
the discussion here will help to stimulate some further experimental and numerical efforts to 
elucidate the relevant vortex dynamics in turbulent flows. In addition, this paper will discuss 
the close analogy between vortex dynamics in classical turbulence and in quantum superfluids. 
A very interesting paper on this subject has been written by Huggins [13], from the vantage 
of the low-temperature physicist. We shall further elaborate upon his discussion, filling in 
some of the details on classical fluid turbulence. Needless to say, comparisons between such 
disparate areas of physics can often be extremely illuminating and can stimulate new questions 
and different points of view in both subjects. 
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2 Vorticity Conservation and Generation at the Boundary 

We begin by briefly reviewing some basic results on vorticity in classical fluid dynamics, its 
associated conservation law and its generation at solid walls. This is necessary in order to set 
the stage for our following discussion on wall-bounded turbulence. Furthermore, even this very 
classical subject is not without controversy, as there is currently some disagreement about the 
proper form of the vorticity source at the wall |15| [T7] . 

2.1 The Helmholtz Equation as Conservation Law 

The evolution of vorticity in classical incompressible fluids is described by the equation of 
Helmholtz [2] 

ata; + (u-V)cj = ((x;-V)u + i/ Aa; + Vxf, (1) 

supplemented by terms representing viscous diffusion with kinematic viscosity v and the curl 
of a non-potential body force f (e.g. from electromagnetic stirring [21] or stress from polymer 
additives |22| [23]). Because this equation is obtained by taking the curl of the momentum 
equation, or incompressible Navier-Stokes equation, it must take the form of a local conservation 
law. Explicitly, ([T| can be rewritten as |14] 

dtu + v-s; = (2) 

or dtujj + diTiij = 0, where 

T,ij = UiUJj - UjUJi + u - + eijkfk- (3) 

This anti-symmetric tensor represents spatial flux of the jth component of vorticity in the ith 
coordinate direction. It plays a fundamental role in all of our following considerations. The 
flrst term on the righthand side of ([3| represents convective transport, while the second term 
represents spatial transport by vortex-shearing, the third term describes viscous diffusion and 
the fourth term gives the induced motion of vorticity due to Magnus effect of the body force. 
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Because vorticity arises as the curl of the velocity, u = V X u, the above conservation law 
must be in some sense "trivial". This is the content of Foppl's Theorem 0, which states 
that the total vorticity integrated over the flow domain O must vanish: f^ud^x = 0. This 
result holds in an infinite domain with sufficiently rapidly decaying vorticity at infinity, or in 
a bounded domain with uniformly moving, solid walls (but not, for example, in a rotating 
container). The standard proof is based on the simple identity u: = V'(a;x), which implies by 
the divergence theorem that J^ud^x = jQ^:x.{uj'n)dA for the outward-pointing normal fi on 
the bounding surface dO,. On a solid wall in uniform motion, cj-n = 0, because of stick boundary 
conditions for the velocity field. There is also no contribution from boundary points at infinity, 
if |a;| = o(|x|~'^) as |x| oo. This result shows that complete cancellation of vorticity must 
occur when integrated over the entire domain of such flows, despite strong localized vortex 
structures of various orientations distributed throughout the fluid. 

2.2 Lighthill's Theory of Vorticity Generation at Solid Walls 

If the vorticity conservation law Q is integrated over the flow domain, one obtains 

^ / u:d^x= [ -EhdA. (4) 

Of course, under the conditions of Foppl's Theorem, which we hereafter assume, Jq^ Xln dA = 0. 
As argued by Lyman [ISj, the relation Q motivates one to consider a = Xln as a source density 
for vorticity at the wall, with strength \cr\ in the direction a = cr/\cr\. As a consequence of 
anti-symmetry of XI, the vortex lines generated are parallel to the surface, since cr-h = |16] . 
However, Wu and Wu [17 have pointed out that such a global argument is insufficient to identify 
cr as a local source of vorticity, since it may be altered by adding any perturbation 6cr such that 
Idn = 0. In the Appendix A, we present a proof based upon the Kelvin Theorem which 

shows that cr = Xlfi does represent a localized source of vorticity at the boundary. This holds in 
the following precise sense: If one considers any loop C in the fluid with part of its length lying 
in the boundary 50, then the contribution of the boundary segment to the time-derivative of 
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the vorticity flux through C is given by 

/ a ■ {nxi)ds 

Jcndn 

where t is the unit tangent vector along the curve C, for a specified orientation. With this 
convention the fiux is measured positive in the direction to the righthand side of the curve C 
as one moves along it, i.e. parallel to nXt. 

Using the expression ^ for it foUows that 

a = -hXvCVxu) +hxf (5) 

due to the no-through condition u-n = for an impermeable wall and the result cj-n = for 
a solid wall in uniform motion. Further insight and simplification can be obtained if one uses 
the Navier-Stokes momentum balance Dtu = — — i/Vxu; + f to give [TH] 

a = nx(Vp + Au). 

This form of the vorticity source density was first derived by Lighthill [24j and Morton |25j 
who obtained the terms from the tangential pressure gradient and the tangential boundary 
acceleration, respectively. The term identified by Morton, nx(Dtu), contributes only if the 
boundary is accelerating continuously or impulsively. If we restrict our attention to constant- 
velocity (or stationary) boundaries, then we recover the original Lighthill (1963) result 

<T = nx(Vp). (6) 

This formula has several remarkable implications [25i. First, one can see that the generation of 
vorticity is essentially inviscid. Second, the generated vortex lines on boundary dO, are parallel 
to the surface isobars, or lines of constant pressure p. 

2.3 Mean Flux of Vorticity in Pressure-Driven Flows 

In turbulent fiow, the vorticity balance ^ may be considered in Reynolds-averaged form: 

atCJ + V-S = (7) 
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where the mean vorticity flux tensor is 



^ - ^ J + ^ijklk- (8) 
The first term UiLOj — UjuJi represents the transport of mean vorticity by the mean flow. On the 



other hand, the term UiUjj — Uj^i represents the transport of mean vorticity by the fluctuations 
and is the analogue of the Reynolds stress in the ensemble-averaged momentum balance. 

We now derive a crucial relationship between mean vorticity transport and mean pressure 
gradients. The Navier-Stokes equation may be written as dtu^ = \(-kirn^im~dk [p + U + 
in terms of S and the potential U of any conservative body force (e.g. gravity). In statistically 
stationary turbulence, the time-derivative averages to zero, so that one obtains 



dk{p + U = ]^ekim^im, ^ij = <^ijkdk{p + U (9) 

To our knowledge, this relation was first derived by Taylor [3^ for a general turbulent shear-flow, 
assuming strict two-dimensionality of the fluctuations. His work is discussed in Section |3]below, 
on turbulence in pipes and channels. More than thirty years later, this result was rediscovered 
by Anderson [lOj, who considered inviscid flow with no body forces. His work will be considered 
in Section |4] on quantum superfluids. The full relationship Q was derived by Huggins |14] . as 
a "Josephson- Anderson relation" for classical turbulence. 

Another useful form of this relation can be derived by using the standard vector-calculus 
identity (u'V)u = uXo; — V(^|up) to write the Reynolds-averaged Navier-Stokes equation as 

— _ 1 — * 

DtUk = -eklmX'im - dkP*, 

with 

% = ^iH^ + - S) + ^--^^ (10) 

and p* = p + U + k, where k = ^ju'P is the turbulent kinetic energy and Dt = dt + u-'V is the 
material time-derivative for the mean flow. Thus, in cases where DfU = 0, 

dkP* = ^^klm^*im^ ^ij = eijkdkP*- (11) 



In this version of the relationship, the vorticity transport by the mean flow is omitted and 
gradients in the "turbulent enthalpy" p* = p + U + k are considered. 

We now consider the implications of these results for turbulent channel and pipe flows. 

3 Channel and Pipe Flows 

As usual, channel flow refers to fluid motion between two parallel, inflnite plane walls separated 
by distance 2h, with a pressure gradient parallel to the planes. Our notations follow those of [1]. 
Thus, we take x to be the streamwise direction (along the pressure gradient), z the spanwise 
direction, and y the wall-normal or cross-stream direction, with y = at the bottom plate. 
The velocity fleld is then represented as u = (u, v, w) in terms of its streamwise, wall-normal, 
and spanwise components. It is assumed that the turbulence is fully developed, in the sense 
that it is statistically stationary and that all averages are independent of z and also (except for 
the mean pressure p) independent of x. The only non- vanishing component of mean velocity is 
u(y), with a maximum value Uc = u{h) at the channel centerplane. The mean vorticity has only 
a spanwise component, iOz{y) = —du/dy, which is negative for y < h and positive for y > h. 

Pipe flow refers instead to fluid motion down an inflnitely long, cylindrical pipe with circular 
cross section of radius p, driven by a pressure gradient down the pipe axis. We use standard 
cylindrical coordinates {r,6,z), so that z represents the streamwise/axial direction, the az- 
imuthal direction, and r the radial/ wall- normal/cross-stream direction. We write the velocity 
as u = {ur, ug, Uz) and similarly for other vector quantities. In fully developed pipe turbulence, 
all average quantities are independent oi t,9 and z, except that mean pressure p decreases 
linearly with z The only non-vanishing component of mean velocity is Uz{r), with a maximum 
value Uc = Uz{0) at the center axis. The mean vorticity has only an azimuthal component, 
^eif) = —duz/dr > 0. Because these two problems are so similar, we shall give full details only 
for channel flow and then comment upon the relevant differences for pipe flow. 
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3.1 Mean Vorticity Flux and Energy Dissipation 



In channel flow, the relations (11) become 



dujz _ y* _ dp* 

dy ~ dx 



V'^'Z - W'Uj'y - U— =T.y,= — (12) 



w'uj'^ - u'uj'^ = S*^ = — (13) 



U'Uj'y - V'iV', = T.,y = 0. (14) 

Note that Ejj = S*^- for all i,j except Tjzx = —^xz = —uljjz + S*^. As should be clear from the 
general derivation, these relations are just the usual momentum balance equations written in a 
different form. For example, it is easy to check that 

d —f—j duJz d'^u 

dy ' dy dy^ 



V'u;'.,-W'u'y = --U'v', -y^ = v (15) 



See also |T], eq.(3.3.13), or [5]. However, now (12)-(14) are seen to have another, equally impor- 



tant interpretation: they relate the pressure gradients in various directions to the corresponding 



cross-fluxes of a transverse component of vorticity. The flrst relation ( 12 ) was already derived 
by Taylor |3] in a slightly modifled form for two-dimensional flow, as mentioned earlier. Taylor 
employed this relation, which is exact, to motivate his well-known "vorticity transfer hypothe- 
sis," which assumes gradient-transport of vorticity. As discussed in [1 , section 3.3, the latter 
hypothesis ignores vortex-stretching interactions and has limited validity (but see j9] for a re- 
flnement). As a consequence, this paper of Taylor's has fallen into a bit of disrepute and it is 
not usually noted that a signiflcant exact result was derived there. 



Of the relations (12)-(14), the most important is indeed (12), since it relates the cross- 
stream flow of spanwise vorticity Yjy^ to the streamwise pressure gradient dp/dx. Note that 
dp* /dx = dp/dx, since U is zero and k is independent of x. It is well-known that the mean 
streamwise pressure-gradient is constant across a cross-section of the channel, i.e. dp/dx = —7* 
independent of y. However, in standard derivations (e.g. [Ij, section 5.2) this appears to be an 
"accident" . Here it is seen that the fundamental reason is the conservation of vorticity: 



It is easy to determine its constant value by the balance of stress at the two walls against the 
pressure-gradient applied over a channel cross-section, giving — 2n^ = (2/i)||, or || = —7=,, = 
—u1/h in terms of the friction velocity u^^. 

The physical picture that emerges is very simple. In view of Lighthill's relation mean 
spanwise vorticity ojz is generated at the two walls, at a rate given by dp/dx = —7*. Due to 
the change in the direction of outward normal vector n, this vorticity is negative at y = and 
positive at y = 2h. The negative vorticity at the bottom wall is transported upward and the 
positive vorticity at the upper wall is transported downward, leading to a constant negative 
mean vorticity flux = —7*. Although the actual geometry and dynamics of vortex lines 

below), what happens "on average" is that straight 
vortex lines are generated in the spanwise direction of equal strengths of vorticity, negative at 
the bottom and positive at the top. These then "cascade" toward the channel centerplane, 
where they annihilate each other. 

It is interesting to compare the above with the more conventional picture of channel flow 
as a "momentum cascade" [T]. Integrating the mean x- momentum balance once in y gives a 
standard formula for the total stress, Reynolds stress plus viscous stress: 

T., = .V-.| = -n^(l-f). (16) 

For < y ^ /i, it follows that Txy ~ —u1 (and likewise Txy ~ u1 for < 2h — y <^ h), leading 
to the notion of a constant momentum flux range. In fact, however, Tiyz = —dyTxy = —u1/h, 
so that the stress can never be strictly constant in any range of y, whereas the vorticity flux is 
exactly constant T,yz = —7* over the whole channel cross-section. 

The importance of the cross-stream flux of spanwise vorticity is that it is directly related to 
energy dissipation. By considering the balance equations of energy in the mean flow (1/2) |mP 
and energy in the fluctuations k = (l/2)|u'P, it is easy to show that 



is quite complex (see subsections 3.2j|3.3 



dy 



(17) 
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Here Ub = ^ Jq^ u{y) dy is the bulk velocity, averaged over a cross-section of the channel, and 



e{y) = z^|Vu'p is the "turbulent dissipation" (per unit mass) in the velocity fluctuations. To 



derive (17) we have used the result that turbulence production balances turbulent dissipation 



when integrated over the channel width: 



2h 



2h 



jdu 
dy 



dy 



2h 



2h 



e{y) dy. 



The relation (17) can be rewritten by multiplication with the mass density q as 



(18) 



where Jf, = qui, is the mean mass flux through the pipe and QS^^^ is the bulk-average energy 
dissipation per unit volume. This relation shows the direct connection of vorticity flux with 
energy dissipation. As we shall see in Section |4) ^ is a special case of a far more general 



relation, valid for both classical and quantum fluids. The practical significance of (12), i.e. 



Tiyz = and (|18|) is that both pressure-drop and energy dissipation can be reduced in turbulent 



flow by decreasing the cross-stream flux of spanwise vorticity. Any mechanism that interrupts 
the "vorticity cascade" will lead to turbulent drag reduction. Body-forces from electromagnetic 
forcing [21] and polymer additives [22l [23] are well-known devices. 



In pipe flow, the relations (11) become 



i^d_^^_^^ ^ dp* 

r dr dz 



u'nto' — u'.ijj'a 



-•Qz 



dp* 
dr 



S,., = 0. 



Similarly as before, Sjj = Sj - for all i,j except E 



The relation (|19|), or Ti^e = is the crucial one. The radial-independence of the streamwise 



Bz 



^ze 



-Uz uje + ^ez ^-^cl 



dz 



(19) 

(20) 
(21) 

dp 
~ dz- 



pressure-gradient follows from vorticity conservation: 

d ( dp^ 



dr \dz 



drT. 



r^rd 



0. 
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The constant value ^ = — 7=,, can be inferred from balancing the stress at the wall with the 
pressure-gradient applied across a cross-section of the pipe, i.e. — 27rr • = vrr^ • ^ so that 
II = — = —2ul/r. The cascade picture "on average" is of vortex rings of positive azimuthal 



vorticity generated at the pipe wall and collapsing radially inward to annihilate at the axis. 
3.2 Scaling Theories and Vortex Structures 

The previous considerations have all been exact. We now develop additional insight using 
theoretical assumptions, involving scaling hypotheses and structural modeling. We consider 
explicitly only channel flow although the discussion applies as well to pipe flow. 

3.2.1 Reynolds-Number Scaling Theories 

We first examine the consequences of high-Reynolds number asymptotics. We shall discuss the 
classical asymptotic matching theory that leads to the logarithmic law of the wall, particularly 
as laid out in the review of Panton [26] . Many of the results that we deduce below from this 
theory can be obtained under more general assumptions, e.g. see |27| [28] . The classical theory 
has, however, the advantages of simplicity and familiarity, so that we shall restrict ourselves 
to it here. As usual, one introduces an inner scaling y+ = u^y/u, u'^ = u/u^ = f{y^,Re^:), 



-u'v' /u1 = g{y~^,Re^,) with Re^, = = The stress relation (16) then becomes 



g{y+,Re,) + 4^{y+,Re,) = l 



(9y+ Re* 

The standard scaling hypothesis of the "law of the wall" is that f{y^,Re*:) ~ foiu^) and 
giy'^^Re*) ~ g^^y^) for Re*. S> 1, and, in that case, go{y~^) +/o(y^) = 1- Matching to an outer 
law leads to the prediction of a "logarithmic layer" with 

/o(y+) ~ -lny+ + B, <7o(y+) ~ 1 ^ 



K Ky~^ 

for /i"*" S> ^ 1. The above results are asymptotically valid only in the inner layer, with 
y <^ h. However, as discussed by Panton [26], these inner laws can be combined with outer laws 
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in a standard technique of "composite expansions" to yield high-Re^: asymptotic formulae that 
are uniformly valid over the whole channel width. In particular, a composite expansion for the 
Reynolds stress is 

~ go{y+) - 1^ (22) 



As shown in [26 , Fig. 40, (22) is a very accurate approximation, even at rather low values of 



i?e=K. This leads to an important consequence. Taking the derivative of —u'v' ~ 1 — — 
one finds [26] 

SO that the peak of the Reynolds stress appears within the logarithmic layer at Up ~ \[^^ 

with maximum value —u'v'^ ~ 1 , ^ . This scaling of the location of the peak Reynolds 

stress has been well- verified empirically pOj |30l EI] . 

There are immediate implications for the turbulent vorticity flux, since, by the first line of 



(15), it is the derivative of the Reynolds stress: ^{—u'v') = v'uj'^ — w'oj'y. It follows that 



w'u>!„ > for w+ < yj 

y « ^24) 



v'io'^ — w'io'y < for y+ > Up ■ 

The condition of constant cross-stream vorticity flux which we derived in section |3] states that 

— — duJz / n 

v'uj'^ — w'uj'.. — V— — = — — < 0. 
" oy n 

The viscous term is negative over the whole channel width, since mean spanwise vorticity is 
negative at y = and increasing monotonically to zero at y = h, so that duJz/dy > for all y. 
However, the turbulent transport term is negative only for y+ > Vp ■ These facts have been noted 
before |27| [28] . We restate them here in order to emphasize that mean cross-stream transport of 
spanwise vorticity is dominated for y+ < by viscous diffusion, which gives the contribution 
of largest magnitude. This result does not contradict any of the assumptions in the derivation of 
the log-layer — and indeed is a consequence of that derivation — but it does contradict a common 
belief that "the overall dynamics of turbulent boundary layers is independent of viscosity" [1] . 
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Further insight can be obtained by considering the constant-flux condition in inner-scahng: 

diot 1 



v'uj', — w'ujI, 



y dy+ Re, ' 
with 



v'Loi - w'uj'y ' = g'{y+,Re,), duJ+Jdy^ = -f"{y^,Re,) 



Right at the wall, g'{0,Re,) = —1/Re, and f"{0,Re,) = 0, so that the vorticity transport 
is entirely due to viscous diffusion in close vicinity of the wall. In the range 1 <^ y^ < yp 
the two functions g' , f" are of order unity, so that the small difference —l/Re-t must be due 
to near-cancellation of almost equal quantities. Thus, the viscous and turbulent transport of 
vorticity are both relatively large and the negative residual value of the cross-stream flux is 
due to a slight dominance of the viscous transport. It is only for y^ > yp that the turbulent 
transport dominates and that one can properly think of an "inertial sublayer" dominated by 
nonlinear interactions. For more discussion of these points, see \n\ I28j. 

One of the consequences of these considerations is that the mean velocity profile u{y) plays 
a critical role in the vorticity cascade over the range y < yp, since the viscous flux equals 
u(Pu{y)/dy'^ . Conversely, the statistics of vorticity production and transport at the boundary 
are a crucial factor in determining the mean profile; for instance, see |32l. This may have 
significant implications on understanding the mechanism of polymer drag-reduction |22| |23] . 



for example. As can be seen from the general expression (10), mean vorticity flux for turbulent 
pipe flow with a polymer additive will consist of three distinct contributions: the nonlinear 
turbulent transport, viscous transport, and direct transport by Magnus effect of polymeric 
forces. Given the known drag-reduction effects of polymers, it follows from our discussion that 
wall-normal vorticity flux must also be drastically reduced by polymer additives. However, this 
may occur a priori in several possible ways. For example, for y < yp either the viscous transport 
may decrease in magnitude or the turbulent transport (which is of opposite sign) may increase 
in magnitude. Further research will be necessary to see which of these possibilities is realized. 
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Figure 1: A horseshoe vortex bundle as originahy conceived by Theodorsen (1952). 
3.2.2 Vortex Structure Models 

Our previous discussion has emphasized the importance of understanding the structure and 
dynamics of vorticity in pipe and channel flow. There have indeed been many proposals to 
explain and model the physics of wall-bounded shear flows in terms of typical vortex struc- 
tures. One of the earlier attempts was that of Theodorsen |33], who suggested that the most 
significant structures are "horseshoe vortices" inclined at 45° to the wall, which form by roll 
up of streamwise vorticity in the viscous sublayer. See Fig. [T] Considerable evidence has been 
found in support of horsehoe vortices (also known as hairpin, A and 0, vortices) from both 
experiment and numerical simulation |34| ESI EHl EZl 133 ES]- However, many other alternative 
proposals have been made. To quote from a review of Fiedler [37 : "Thus, when studying 
the literature on boundary layers, one is soon lost in a zoo of structures, e.g. horseshoe- and 
hairpin-eddies, pancake- and surfboard-eddies, typical eddies, vortex rings, mushroom-eddies, 
arrowhead-eddies, etc It is fair to say that no clear consensus has emerged as to which 
vortex structures are the most relevant to the dynamics and statistics of the turbulence. 

Most of the effort in identifying coherent-vortices and in constructing vortex models of the 
turbulent boundary layer has gone into explaining their contribution to the Reynolds stress. The 
importance of maintaining a (nearly) constant Reynolds stress across the "inertial sublayer" has 
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been widely recognized. However, almost no effort has been made to explain the cross-stream 
transport of spanwise vorticity in terms of typical vortex structures. (But see O El El El [9] , 
discussed more below.) Furthermore, the fundamental constraint of maintaining a constant 
mean vorticity flux across the pipe or channel has been almost universally ignored. 

Here we shall discuss such issues in the context of the "hairpin vortex" model of Perry, 
Chong and collaborators |40l \^T\ H2l UHl Ej. This model is one of the most fully articulated 
attempts to explain, in analytical terms, the features of the turbulent boundary-layer by means 
of representative vortex-structures. The Perry-Chong (PC) theory combines the Theodorsen 
proposal [33^ of "horseshoe" or "hairpin" vortices with Townsend's "attached eddy hypothesis" 
|45j . In this picture, the boundary layer consists of a scale-invariant hierarchy of hairpin-like 
structures, which are attached to the wall in the sense that their characteristic lengths are 
proportional to the distance at which the eddy extends above the wall. By averaging over an 
ensemble of such structures, the PC theory yields important quantitative results, such as a 
logarithmic mean veocity profile and constant Reynolds stress in the overlap region, profiles of 
normal Reynolds stresses, tensor energy spectra as functions of distance from the wall, etc. 

PC theory is sufficiently developed that the turbulent contribution to vorticity flux S*"^* = 
v'oj'^ — w'to'y may be calculated by evaluating each term separately, e.g. as integrals in wavenum- 
ber over the model energy spectra. However, for our purposes a simpler approach suffices, which 
is to use again the identity ^{—u'v') = v'uj'^ — w'uj'y. We may start from the PC formula for 
the Reynolds stress as an average over hierarchy lengths 5, e.g. equation (27) in |32]: 



Here 5i = 0{u/u^) is the length-scale of the smallest hierarchy, whereas A^; = 0{h) is the 
length-scale of the largest. The function Iij{y*) is obtained by integrating over spectra of 
the representative vortex structures. See also Townsend [IS], who referred to Iij{y*) as an 
"eddy-intensity function." Both PC and Townsend required a distribution of hierarchy lengths 
Ph{S) = M/5 for 6i < 5 < Ae in order to reproduce a logarithmic mean velocity profile and 




(25) 
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a constant Reynolds stress. We may follow their arguments to study also the derivative of the 
stress, e.g. see |l0], Appendix B. Thus, changing integration variable to y* = y/S 



u'y. = Mul / li.iy*) ^ 

^ Jy/si y 



and then differentiating gives 

d Mul 



E 



(26) 



dy y 

Using the result Ixyiu*) for y* > I and Ixy{y*) ~ -Qy* for y* <.! [10] gives 

— — d,^^^ MQul 

v'ljj — w uo' = —i—u'v ) 

^ dy Ae 

for 5i ^ y <^ Ae- This same result may also be obtained by differentiating the final result for 
the Reynolds stress in Appendix B of |l0]. Since A^; = 0{h), the PC theory yields a constant 



of the correct order, v'u)'^ — w'uj'y = 0(l/i?e*). 

These results appear satisfactory, but a little thought shows that the PC theory is adequate 
to explain vorticity flux only over the range y^ < y^ <^ h~^, in the notations of the previous 
subsection. It is only over that range that the turbulent vorticity flux is constant and, otherwise, 
the viscous diffusion is dominant. To obtain the correct negative sign of the turbulent flux over 
the range yp < y^ <^ (and also to obtain the correct sign of the Reynolds stress [30]), the 
constant Q in the above calculation must be positive. However, as we have seen in the previous 
subsection, the turbulent vortex flux must reverse sign and become positive for < y^ < y^. 
The PC theory, as it is presently constituted, is inadequate to explain this sign-reversal. This 
is not a failure just in the buffer layer but over half of the range of the traditional log-layer. 

Let us consider briefly one possible remedy. As mentioned above, it has usually been 
assumed in the PC model that Ixy{y*) = for y* ^ 1. For example, see Figure 14 of |41j 
or Figure 4 of |32], where Ixy{y*) = for y* > 1. This corresponds to the severe assumption 
that vortex hierarchy elements of length-scale 6 make no contribution to Reynolds stress at 
wall-normal distances y > 5.11 one hypothesizes instead a slow power-law decay of the precise 
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form Ixy{y*) ~ ~y^' y* ^ ^' ■^ith P another universal, positive constant, then (26) yields 



dy \ y2 A 



for (5i ^ y ^ A^;, in qualitative agreement with eq.(23). The agreement becomes exact for a 
suitable choice of the constants M, P, Q. An interesting problem for future investigation might 
be to determine which "representative" vortex structures, if any, lead to the necessary —1 
power-law decay of the eddy-intensity function Ixy{y*) for y* ^ 1. 

3.3 Experiments and Simulations 

It would be very useful to get some additional clues from experiment and simulation. There 
seem, however, to be very few empirical studies of the turbulent vorticity flux or its relevant 



constituent averages, v'uj'^ and w'oj'y. We know of no experimental measurements at all for 
turbulent flow in pipes or ducts. The most exhaustive investigations have been carried out in 
zero- pressure turbulent boundary layers, by Klewicki and collaborators. Those studies include 
wind-tunnel boundary layers with Ree = 1010 ~ 4850 [5l [6] and atmospheric boundary layers 
with Reo = 2 ^ Ax 10^ j7j. We also know of only two studies of these velocity- vorticity corre- 
lations by direct numerical simulation of channel flow, the iJe* = 250 calculation of Leighton 
& Handler reported by Bernard [9 and a -Re* = 200 DNS of Crawford & Karniadakis |i21| . 
The picture which emerges from the lower Reynolds-number studies [5l El [9l [21] is as follows: 



v'uj'^ > > w'oj'y for y+ < 12 



> v'uj'^ > w'oj'y for 12 < y+ < 



> w'u' > v'uj'^ for y > yj, 



A question mark is included in the last line since the DNS of O [2T] indicates that w'tUy becomes 
slightly positive for ^ 50 > ~ 30. The results for y~^ > agree rather well with the 
physical picture of vortex dynamics in the PC model jlQ]: "the A-vortex stretches under its 
own mutual induction with its image to give, in the large, plane strain ... The plane strain 
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brings the legs of the vortex together" . The hfting of the negative streamwise vorticity in the 
head and legs of the hairpin produces v'uj'^ < 0. The simultaneous collapse of the hairpin legs 
towards each other, carrying opposite signs of wall-normal vorticity, produces w'uj'y < 0. The 
slight positive value of w'tOy for larger y values has been interpreted by Bernard [9j also in terms 
of hairpin structures, as due to the reversal of w' at some height above the counter-rotating legs 
of the vortex. Our calculation in the previous subsection indicates that w'uj'y > v'uj'^, yielding a 
net negative vorticity flux. For more detailed discussion of these correlations in terms of typical 
vortex structures, see |6], Section III and |9], Section 4. 

On the other hand, the PC model does not seem to accord well with the lower Reynolds 
observations for < y^. Klewicki et al. [6] have made a detailed quadrant analysis of the 
contributions to v'uj'^ at y^ = 5.3, 14.2 and 26.3 in the BL at Reg = 1010 with y^ =40. 
Although v'lv'^ < for y+ = 14.2 and 26.3, as the hairpin picture might suggest, [B] have found 
that the dominant events are 2nd-quadrant, i.e. v' < and uj'^ > 0. Furthermore there is 
strong dependence on Reynolds number, with the magnitude of v'uj'^ decreasing sharply with 
increasing Reynolds number for Reg = 1010 ~ 4850. 

These tendencies persist in the high Reg data from the ABL in [7j. In fact, their smooth- 
wall data at Reg = 2 x 10^ show v'uj'^ > for the whole range of y'^. Furthermore, w'iUy 
appears to change sign, being negative for y+ < y+ and positive for y+ > y+ (using rough-wall 
data). The sign for y+ > y^ may be due to growth of hairpin vortices in the spanwise direction 
\39\ H6l ITT] . Clearly, further experimental and numerical studies would be most helpful, in 
order to verify the limited data and to map out all the dependencies on wall-distance, Reynolds 
number, roughness, etc. Particularly illuminating would be conditional sampling studies, with 
VISA (variable- interval space-averaging) sampling technique and quadrant conditions like those 
applied to Reynolds stress j351l36j. but now carried out for the vorticity flux. Such investigations 
could help to reveal the vortex structures and dynamics which contribute most significantly to 
the cross-stream transport of spanwise vorticity. 
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4 Comparison of Classical and Quantum Fluids 

We shall now discuss briefly some of the corresponding problems in the field of quantum super- 
fluids. A complete review of superfiuidity and superconductivity would be out of place here and 
we refer the reader interested in further background to standard texts |48| 149] . We just remark 
on a few basic facts, the most important of which is quantization. As first noted by Onsager 
|50j and Feynman [51], the vorticity in neutral superfluids resides entirely in thin filaments 
whose circulation is quantized in units of k = /i/m, where h is Planck's constant and m is the 
atomic mass of the condensed Bose particle. For liquid ^He, k ~ 9.97 x 10"^ cm^/sec and the 
core radius of the vortex lines is of order of angstroms. An analogous fact holds for charged 
superfluids or superconductors, in which the condensed Bose "particle" is generally a Cooper 
pair of two electrons of charge 2e. As first pointed out by London [52j, the magnetic fiux in 
a superconductor is likewise quantized in units of = h/2e = 2.07 x 10~^^ webers. London 
considered what are now called "Type-I superconductors" which are perfect diamagnets and 
expel all magnetic fields (the Meissner effect). Thus, he had in mind the magnetic fiux through 
superconducting rings and other non-simply-connected samples. However, it was pointed out 
later by Abrikosov |53| that there is another class of "Type-II superconductors" which are 
penetrated by magnetic fiux tubes with fiux quantized in units of ^q. 

We now discuss the "phase-slip" phenomenon, which is closely analogous to the results for 



classical turbulence presented in section 2.3 These ideas were first employed by Josephson for 
superconductors and extended further for both superfiuids and superconductors by Anderson, 
hence often called the "Josephson- Anderson relation" |T3]. The classic presentation is in the 
review article by Anderson [101. He argued there that the drop in chemical potential (per 
unit mass) between two points Pi and P2 in a neutral superfiuid should be related to the rate 
dn/dt at which quantized vortex-lines cross an oriented line between the two points: 



/i(P2) - /^(-Pi) = K.dn/dt (27) 
Here the overbar indicates a time-average. In an appendix — entitled "A 'New' Corollary in 
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Classical Hydrodynamics?" — Anderson also pointed out that this relation for quantum fluids 
has a precise classical analogue, of the form: 



This is just the relation (|9| that we derived earlier, in the special case of an inviscid fluid, with 
no body-forces, and in line-integral form. It is a generalization of the relation noted by Taylor 
|3] for shear turbulence. A similar relation can be obtained for type-II superconductors 



now relating the difference in voltage between two points to the rate at which a line from Pi to 
P2 is crossed by quantized flux |55| [56]. This relation was the basis of Anderson's "flux-creep 
theory" to explain resistive behavior in type-II superconductors above a critical current. 

These ideas were perhaps most carefully elaborated by Huggins [11] . His paper is especially 
interesting for us, because he employed the classical fluid equations as his model, despite in- 
tended applications to quantum superfluids. This is legitimate, because a superfluid will obey 
the laws of classical hydrodynamics except close to the cores of vortex lines. Huggins makes 
several assumptions that are adapted to the superfluid problem, e.g. that all the vorticity is 
conflned to compact vortex tubes with total circulation of magnitude k. However, almost all of 
his discussion can be easily generalized to the case of continuous vorticity distributions of any 
magnitude (see Appendix [B]) . 

The starting point of Huggin's treatment is Kelvin's minimum kinetic energy theorem |57] . 
As is well-known, Kelvin showed that, among all incompressible velocity flelds with the same 
boundary conditions in a flnite domain $7, the one with minimum energy is the potential flow. 
To prove this, Kelvin decomposed the velocity fleld as u = -|- where u,^ = V</> is the 
potential flow with the specifled inflow-outflow boundary conditions and u^^ is a fleld containing 
all of the vorticity of the fluid but with no flow through the boundary. In that case the kinetic 
energy decomposes into a sum of kinetic energies for the potential and vortex flelds, so that 




(28) 




(29) 



21 



E = + > E^. Huggins discusses a classical fluid model in the form 

a^u = -V/i + uXa; + g, (30) 

where the enthalpy h = p+[/+^|up includes the contribution of any potential energy U whereas 
g contains all the non-potential forces, e.g. g = — z>(Vxa;) + f if there is only a viscous force 
and a non-conservative body-force. Huggins further considers flow through a channel driven by 
a difference in the enthalpy at the inflow section and enthalpy /loni at the outflow section. 
In this setting [TT] proves that 



E = J{hin-hout)+ / QVi-^d^x (31) 
E^ = J{hin — hout) — i^J(f, (32) 

E'o; = kJ,^ -I- / QVL-^(fx (33) 

Jo. 

Let us consider each of these relations in turn. 

The first result ( |31| ) is just the standard energy balance relation. Since J is the total mass 
flux through the pipe, J{hin — hout) equals the energy input into the fluid through work by 
the applied enthalpy difference. On the other hand, D = — i?U'g d'^x just represents the 



dissipation by the non-potential forces. The second relation (32) is more interesting. This 
result shows that the energy supplied by the enthalpy difference goes entirely into the potential 
flow. What acts as "dissipation" on the potential flow is the term kJ^, which transfers energy 



to the vortex system. As seen from the final relation (33), this energy is then dissipated by the 
non-potential forces, e.g. ordinary viscosity, acting on the vortices. The transfer term is the 
product of the quantum of circulation n and the rate at which vortex lines are crossing the 
mass flux in the potential flow. For example, for a single length of vortex line 



K,Jff) = K I QVi(j,-VivXd\ (34) 

J vort 

The quantity u„ is the velocity of the element d\ of vortex line, so that VL^xdX represents the 



rate at which the line-element sweeps out area. It is easy to recover Anderson's relation (27) [or 
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(28)] from (32). For example, if vortex lines are crossing the entire mass flux J in the channel 



at a frequency then Js = Jv. In that case, hin — hout = i^^^ for steady-state flow. 



Note that it is crucial in (34) that u„ ^ u^, or otherwise J^f, = 0. Thus, it is important that 
quantized vortex lines are not material lines, moving with the local potential velocity u^. An 
analogous fact is true for vortex lines in classical flow which, because of the non-potential force 
g, do not move with the local fluid velocity u. In general one may decompose the non-potential 
force into components transverse and longitudinal to the local vorticity vector: 

% = /\u.Xu + au, (35) 

where Au = ujXg,/\u}\'^ , a = a;'g/|a;p. The vector Au can be interpreted as a "drift velocity" 
of vortex lines with respect to the local fluid velocity, so that the vortex line moves with net 
velocity = u -|- Au. In that case, the transverse force has the form of a "Magnus force" 
(u^ — u)xa;, similar to that which acts on superfluid vortices |58| I59j . On the other hand, 
there are non-vanishing longitudinal forces as well in classical fluids and these lead to non- 
trivial effects, such as helicity cascade [60^ and kinematic a-effect [61 . The effect of both terms 
on vortex motion is clearly seen in turbulent flows, where viscosity leads to quite different 
properties of material lines and vortex lines |62J- These effects appear to persist even in the 
zero- viscosity limit, when vortex lines behave as material lines only in an average sense |63| 164] . 

It is interesting to note that the phase-slips of quantized flux lines and their consequent 
dissipation lead to a "drag reduction problem" in superconductor technology. In type-II su- 
perconductors the result analogous to Huggins' above is that crossing of supercurrent by flux 
lines at a rate produces an energy dissipation I{Vi — V2) = ^olif), where / is the total elec- 
tric current. This motion of vortices produces a "quasi-resistance" in type-II superconductors 
above some critical current Ic, which destroys their useful superconducting properties. This 
problem is particularly severe for the newer class of high-temperature superconductors (HTS), 
since the vortices there are subject to weaker restoring forces and stronger thermal fluctuations. 
Unlike the vortices in conventional type-II superconductors which are frozen into an Abrikosov 
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lattice or vortex crystal, the vortices in HTS melt at high enough temperatures into a mobile 
"vortex liquid" that produces sizable resistivity. See [651 EHl E7] and, for an excellent popular 
account, [68.^. One solution to this problem is to "pin" the vortices to prevent their motion. 
Methods that have proved successful in more conventional type-II superconductors are to dope 
the superconductor with impurities or to manufacture wires and films from sintered powders 
(e.g. see |69l [70] as representative papers). Anything that prevents free cross-current motion 
of the quantized vortices restores the superconducting properties and permits electric current 
to flow with no resistance. 

An even closer analogy exists between classical turbulence and turbulence in neutral super- 
fluids, as originally conceived by Feynman [51]. To quote his own words: 

"In ordinary fluids flowing rapidly and with very low viscosity the phenomenon of 
turbulence sets in. A motion involving vorticity is unstable. The vortex lines twist 
about in an ever more complex fashion, increasing their length at the expense of 
the kinetic energy of the main stream. That is, if a liquid is flowing at a uniform 
velocity and a vortex line is started somewhere upstream, this line is twisted into 
a long complex tangle further downstream. To the uniform velocity is added a 
complex irregular velocity field. The energy for this is supplied by pressure head. 
We may imagine that similar things happen in the helium." — Feynman (1955) 

Experiments on turbulent flow of superfluid ^He in circular pipes [71] and rectangular ducts |72) 
show similar behaviors as classical turbulence, both in drag resistence as functions of Reynolds 
number and in mean velocity profiles. In addition to the complex, random motion of vortex 
lines discussed by Feynman above, there must also be organized motion. Just as in the classical 
case, there must be a "cascade" of azimuthal vorticity along the radial direction toward the pipe 
axis. As follows from our discussion above, this is essential for the vortex lines to extract energy 
from the pressure head. Vortex structures like "hairpins," "rings," etc. will also doubtless form 
and contribute to the necessary cross-stream flux of vorticity. 
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5 Discussion and Conclusions 



We have presented in this paper a systematic picture of pressure-driven turbulence in pipes and 
channels as spatial "inverse cascades" of spanwise vorticity in the cross-stream direction. This 
complements the more conventional picture of these flows as "forward cascades" of streamwise 
momentum, also across the stream but toward the viscous layer at the wall. Following Huggins 
|14j . we have shown that the spatial flux of vorticity is necessary to produce a pressure-drop 
down the pipe and consequent dissipation. This flux is dominated by viscous diffusion of 
vorticity at distances to the wall less than the peak of the Reynolds stress, well into the classical 
log-layer. The Perry-Chong model based on "representative" hairpin/horsehoe vortices has been 
shown to be inadequate to describe the cross-stream vorticity flux over that range, but viable 
at distances greater than the Reynolds-stress maximum. We have shown that similar results 
hold for quantum superconductors and superfluids, and, especially, for superfluid turbulence 
driven by applied pressure gradients in pipes and channels. 

We hope that one consequence of this work will be to stimulate a greater effort by ex- 
perimentalists and numericists to discover the detailed vortex structures and dynamics which 
contribute to the cross-stream vorticity transport, as a function of distance to the wall and of 
Reynolds number. Theorists should also recognize the condition of constant wall-normal vor- 
ticity flux as an important constraint on statistical models and closures. On the practical side, 
we cannot emphasize too strongly that pressure-drop and cross-stream transport of spanwise 
vorticity are mathematically equivalent. This is a very universal fluid-dynamical result, appli- 



cable not only to turbulence but also to laminar and transitional flow. See eqs. (32), (33) or the 



very general formulation in eq. (39) of Appendix B, for instance. In fully-developed turbulence 
in channels and pipes, which enjoys all the statistical symmetries permitted by the boundary 
conditions, this dynamical constraint leads to constant time-average flux of spanwise vorticity 
across the main flow. The laminar solution in channel geometry with parabolic velocity profile 
has also a constant vorticity- flux, given entirely by viscous transfer (see |14j, section 13). In 
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the turbulent case, there is an additional mean contribution from nonlinear transport. For a 
fixed mass-flux, the vorticity transfer in turbulent flow is much greater than in the laminar 
state, leading to a greatly enhanced streamwise pressure drop and enhanced energy dissipation. 
Any method for reducing drag must therefore reduce the cross-stream transfer of vorticity. It is 
hoped that the vortex-cascade framework presented here will provide new perspectives on the 
difficult problems of wall-bounded turbulence. 

A fresh view of the role of the vorticity may be of use, in particular, in the problem of 
polymer drag reduction. It has long been suspected that the critical action of polymers in 
turbulent drag reduction is through their interaction with vortices |22l l23l [73] . The same is also 
true for various other drag-control strategies, such as with electromagnetic body forces [21]. A 
persistent puzzle has been to understand how polymer-induced changes in vortex structure are 
mediated into substantial changes of the large-scale structure of the flow. The vorticity flux 
tensor (|8| is a natural locus for critical polymer activity. There is generally a scale separation 
in high Reynolds number turbulence between the peaks of the spectra of the velocity and the 
vorticity. However, jTj have found that the velocity-vorticity cospectra in boundary layers are 
not generally dominated by intermediate wavenumbers between these two peaks, but instead 
get dominant contributions from wavenumbers near the peak in the associated vorticity spectra 
and, in some cases, also near the peak in the associated velocity spectra. Thus the velocity- 
vorticity correlations that contribute to the turbulent vorticity flux are a natural bridge between 
the small-scale and large-scale features of the flow. 
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A Appendix: The Vorticity Source Density at a Wall 

It is the purpose of this appendix to derive the fohowing expression 

(T = n' X [uV X u — f] 

for the vorticity source density at the wall. Here, in contrast to the main text, we have preferred 
to use the inward-pointing unit normal h' = —h at the boundary dO, of the domain. Let us 
consider any unit vector e-j locally tangent to the surface, n • ej = 0, and consider the small 
square loop C normal to e^, with one edge lying in the boundary surface dO,: 




Figure 2: A test loop orthogonal to the wall. 

Here, x ej = = n' and the edge tangent to e^- lies in the surface d^l. If S is the square 
surface bounded by C, 

/ uJidxjdxk = f u ■ dx. 
Js Jc 

Letting S{t) and C{t) be the surface S and loop C materially advected by the fluid, then by 
the Kelvin Theorem 



dt 



f uJidxjdxk\t=o = 4; / • (ix|t=o = / [-z^V x c^; + f] • dx. 

JS{t) JC{t) JC 
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The bottom segment parallel to — which is stationary in the surface because of the stick 
conditions — contributes a term 

/ [-z^V X cj + fjjdxj 
Jcndn 

to the time-derivative of the flux JgUJidA. We therefore identify 

CTj = -i/(V xu:)j + fj 

or, for any direction, 

a = n X [uV X — f], 

in agreement with our earlier finding. 

The above argument also reveals the precise meaning of the "vorticity source density" cr. If 
one considers any loop C with part of its length in the boundary dQ, then the contribution of 
the boundary segment to the time-derivative of the vorticity flux through C is given by 

/ <T • (t X h')ds 

Jcndn 

where t is the unit tangent vector along the curve C, for a specified orientation. With this 
convention the flux is measured positive in the direction to the righthand side as one traverses 
the curve C, i.e. parallel to t x n'. 

The expression that we have derived for the "vorticity source density" 

CT = n' X [z^V X — f] 

is the same as that obtained by Lyman [15]. There seems, however, to be some controversy 
about this result. It has been argued by Wu & Wu [T7] that Lyman's prescription (and ours!) 
is incorrect and they suggest that fi' x z^(V x u) be replaced instead with 

'^Wu = = • '^)^ 

with outward-pointing normal n. The two expressions have equal integrals over the boundary 
dQ of the flow-domain since, by the divergence theorem and vector-calculus identities, they 
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both equal uV'^uj d?x. The two differ locally by a term ( Vn)a;, which vanishes for a flat 
surface. In particular, the two expressions agree for channel flow. However, with the precise 
interpretation of the vorticity source density given by us above, our result is a direct consequence 
of the Kelvin Theorem. We note furthermore that a ought to lie parallel to vorticity vectors 
generated within the boundary surface and thus satisfy fi-cr = 0. While Lyman's prescription 
always satisfies this condition, that of Wu & Wu does not in general. 

B Appendix: The Detailed Josephson Relation 

We present here a proof of the "detailed Josephson relation," due to Huggins Although all 
of the results presented below are essentially contained in his paper, the generality of those re- 
sults is somewhat obscured by various special assumptions invoked in that work. Our derivation 
has a slight novelty and hopefully makes the main ideas more clear. 

We consider a generalized channel geometry, where the simply-connected flow domain $7 
has boundary consisting of an in-flow surface Sin , an out-flow surface Sout , and wall surface S^ ■ 
The boundary conditions on the velocity u are: 

= = ^out, u\s^ = 0. (36) 

Huggin's argument builds upon the proof of the Kelvin minimum-energy theorem [57]. Thus, 
the velocity is decomposed as u = -|- u^^, where u,^ = V(/> is the unique potential flow field 
satisfying the incompressibility constraint V-u^ = = and the through-flow boundary 
conditions 

n-u^L, = n-Ui„, n-u^L = ii-Uout, n-u<AL =0, (37) 

inherited from the full velocity u. Thus, u^j = u — u,^ is also incompressible, carries the full 
vorticity of the flow = Vxu^^), and satisfies the no-through b.c. n-u^lg^ = 0. It follows 
from these facts that 

Jo. Jn Jdn 
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This orthogonahty of and is, of course, the essence of Kelvin's theorem. 



One can now substitute = ^(f) into the equation of motion (30) to obtain an equation 
for the time-derivative of : 

Uo; = -V/i' + uxo; + g 

with h' = h + (j) = p + U + ^|up + (j). This result may integrated along any path C in the fluid 
from point Pi to P2, yielding 

rP2 rP'2 

h'{P2)-h'{Pi)+ u<^.(ix= / (uxcj+g)-(ix. (38) 
JPi JPi 

This result is already close to the desired relation. For example, if one Reynolds-averages with 



respect to a stationary ensemble, then u^^ = and thus (38) reduces to Anderson's relation 



(28), in a slightly generalized form. 

The key observation of Huggins is that another kind of average, a "mass-flux average," also 
makes the time-derivative term vanish. To explain this average, we must say a few words about 
"intrinsic" or "natural" coordinates for the incompressible potential flow u^. This system takes 
as coordinate surfaces the iso-potential surfaces 5^ and the streamlines of u^. As a coordinate 
variable along streamlines one may take (j) itself or alternatively arclength related by dcf) = 
u^M. The line- vector element along streamlines is thus M = u,^ di. On a selected iso-potential 
surface one may introduce curvilinear coordinates to label the individual streamlines, 
completing the set {(f), ip, x) of natural/intrinsic coordinates. The area element on is given by 
dA = I X ^\dtpdx, in terms of which another positive measure may be defined as dJ = gu^f^dA. 
This measure represents the mass-flux carried by the potential flow u,^ across each subset of 
the surface 5*0. Note that it does not matter in the definition of dJ which iso-potential surface 
is considered, since dJ is constant along streamlines by incompressibility. For a general surface 
transverse to the potential flow, not necessarily an equi-potential surface, dJ = g\U(f,-h\dA. The 
entire construction simplifies in two space dimensions. In 2D intrinsic orthogonal coordinates 
are (0, "0)) where ip is the streamfunction which is the harmonic conjugate to the potential (p. 
In that case dJ = gdip. 
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The derivation is now easily completed. Take the integral in (38) along the streamlines 



crossing from surface Sin to surface Sout, and average with respect to dJ. It then follows from 



the b.c. (37) that the mass-flux of the potential flow is equal to the mass- flux of the full flow. 



The key term to analyze is / dJ f u^j-di. Using gU(f,d^x = QU(f,dA • u^pdi = dJ d£, this may be 
rewritten as 

dJ / Uuj-di = Uuj-UA)d^x = 0. 



The latter integral vanishes by the same reasoning used to prove orthogonality of u,^ and u^. 
We finally obtain the detailed Josephson relation of [H] : 



J dJK^-J dJK^t = -l dJ Jiuxu + gyde 



(39) 



The line-integrals on the righthand side are taken along streamlines of the potential flow and 



represent the generation of circulation along those lines. Substituting the decomposition (35) 



for g one can see that there are two essential contributions, one from the line-integral of u^Xu: 
which represents the transverse motion of vorticity across streamlines at velocity u^, = u -|- Au 
and the other from the line-integral of au which represents acceleration by longitudinal forces. 
Note that one may also write 

uxu + g)-de = ^ / eijkT,ijdek, 



in terms of the vorticity flux tensor (|3| . The lefthand side of ( 39 ) may be rewritten as 

J dJ /i-„ - j dJ Kut = J [{h'ijj - {Kut)j] , 
where J is the total mass flux and = (1/J) J{-)dJ defines the mass-flux average. The 



energy relations (31)-(33) may be similarly generalized, with J{hin — hout) in (31) and (32) 



replaced with J[{h, 



in I J 



{hout) j] , and the transfer term kJ^ in (32) and (33) replaced with 



the righthand side of (39). For details of the proofs, see jllj . 
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